Magnetic and Quasiparticle Spectra of an itinerant J\ — J 2 Model for Iron Pnictides 
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We calculate the magnetic and quasiparticle excitation spectra of an itinerant Ji — J2 model for 
iron pnictides. In addition to an acoustic spin-wave branch, the magnetic spectrum has a second, 
optical branch, resulting from the coupled four-sublattice magnetic structure. The spin-wave velocity 
has also a planar directional anisotropy, due to the collinear/striped antiferromagnetism. Within the 
magnetically ordered phase, the quasiparticle spectrum is composed of two Dirac cones, resulting 
from the folding of the magnetic Brillouin zone. We discuss the relevance of our findings to the 
understanding of both neutron scattering and photoemission spectroscopy results for SrFe2As2. 
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Since the discovery of superconductivity above 50 K in 
RFeAsOi-zF^ (R = La, Ce, Sm, etc.) the term high 
temperature superconductivity can no longer be associ- 
ated exclusively to cuprates. The new iron-based pnictide 
compounds, including the MFe2As2 (M = Ba, Sr, etc.) 
family exhibit, just like cuprates, a layered antiferro- 
magnetic (AF) structure in the parent compound, which 
gives room to superconductivity upon doping [3j. One 
important difference, however, is related to the transport 
and optical properties of their parent compounds. While 
cuprates are Mott antiferromagnetic insulators, pnictides 
are metals, albeit not very good ones, with a very rich 
multi-band Fermi surface [4j. Despite this difference, 
it is believed that, in both cases, strong correlations of 
strength U among transition metal electrons, from Cu in 
the case of cuprates and from Fe in the case of pnictides, 
play a very important role jf|. In fact, it has been ar- 
gued that, at zero temperature, T = 0, these two systems 
could be represented as very closely, but at oposite sides, 
to a Mott metal insulator transition, in the T x U phase 
diagram This would explain, for example, the experi- 
mentally observed strong renormalization of the coherent 
quasiparticle spectral weight at the Fermi level in pnic- 
tides [6j , which would be shifted to form incoherent lower 
and upper Hubbard bands, just like in prototypical Mott 
insulators such as cuprates. 

The incoherent part of the spectrum in iron pnictides 
can be well described in terms of localized moments at 
the Fe positions, interacting via a superexchange J\ be- 
tween nearest Fe neighbors, Fe— Fe, and a second su- 
perexchange J2 between As bridged next-to-nearest Fe 
neighbors, Fe— As— Fe [7]. Alternatively, an itinerant 
description, in terms of coherent multi-band electrons 
and holes, Fermi surface nesting, and spin-density-wave 
(SDW) instability, can also be used for the understanding 
of both the antiferromagnetic and superconducting states 
in iron pnictides [1, 0|. In this letter, we shall adopt a 
combined description, merging these two degrees of free- 
dom, in which local magnetic moments are coupled, via 
Hund's exchange interaction [Io| . to itinerant quasiparti- 
cles, to be described by a minimal two-band model ll|. 



Let us start by analyzing the magnetic excitation spec- 
trum. One important difference between having S = 1/2, 
as in cuprates, or S = 1 (or higher), as in pnictides, is 
the possibility of single ion anisotropy, in the later case, 
which arises from relativistic corrections and the spin or- 
bit coupling. The simplest model that captures all the 
essential local moment physics described in the previous 
paragraphs is the extended J\ — J2 spin-Hamiltonian 

h = j x ^ ■ sj + j 2 s < • § i - K Y,&) 2 ' c 1 ) 



where 3\ > and J2 > are, respectively, the antifer- 
romagnetic superexchanges between nearest-neighbors, 
and next-to-nearest neighbors, spins on 

a two dimensional square lattice, and K > is the single 
ion anisotropy coupling constant. 

The existence of two superexchanges J\ and J2 in 
the spin Hamiltonian ([T]) renders the antiferromagnetism 
collinear, with wave vectors at (n,0) and/or (0,7r) [T^ |. 
This is in agreement with inelastic neutron scattering in 
(Sr,Ba)Fe 2 As 2 [H, [3, which exhibit peaks at the AF 
zone center Q = (1, 0, 1), for energies between 5—15 meV. 
This situation is markedly different from a Neel ordered 
state with wave vector at (n, 7r), as found in cuprates, 
where a single superexchange J between Cu— O— Cu is 
present. As discussed in the literature, such wave vector 
degeneracy in the classical ground state of the Hamilto- 
nian ([T]) gives rise to an extra Ising symmetry (l5- 17|. 
which is broken at a different temperature than the AF 
ordering one. The hierarchy of symmetry breaking as 
T is lowered is: paramagnetic (PM) /Ising symmetric, 
PM/Ising broken, and AF/Ising broken phases (l6| . 

Previous linear spin wave studies of the J\ — J2 model 
with K = have considered, as a starting point, a classi- 
cal ground state with a two sublattice structure, and have 
obtained that the quantized magnetic spectrum would be 
composed by a single gapless spin- wave branch [l8[ , with 
7kj(k) ~ |k| for small wave vector k. The existence of 
such gapless mode reflects the full rotational invariance 
of the spin Hamiltonian (JIJ, when K — 0, and is in agree- 
ment with Goldstone's theorem. We shall refer to this 
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branch as the acoustic branch. On the other hand, inelas- 
tic neutron scattering revealed the gapped nature of the 
magnetic spectrum [13l Il4| , indicating that the single- ion 
anisotropy parameter K ^ 0. However, the absence of a 
clean step in the neutron integrated intensity (expected 
for the case of a single gap), makes room for the exis- 
tence of two or more magnetic branches IH 14 1 . As we 
shall soon see, the actual jour sublattice structure of the 
AF/Ising broken classical ground state, with two cou- 
pled, interpenetrating Neel ordered states, gives rise to 
a second, optical spin- wave branch (already for K = 0). 
Below we will consider K ^ 0, when rotational invariance 
is broken and all branches become gapped. 

In order to understand the origin of the optical mode 
in the J\ — J2 Hamiltonian fl} it suffices to consider the 
isotropic case, K = 0. For J2 ~ 2Ji, as the experi- 
ments suggest, the ground state is a collinear antiferro- 
magnet, composed by two coupled, interpenetrating Neel 
ordered states We span spins in a coherent basis for each 
of these Neel states, which we label as A and B, and we 



write S? = SSI? = S 



and Sf = Stl? = S 
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where Q = [it / a, it / a) is the ordering wave vector for 
each Neel state, t\a,b and La,b are, respectively, the 
staggered and uniform components of the spins belong- 
ing to the two states, and s — S/a d is the density of spin 
in the unit cell. We use the constraint B = 1. 

After integrating out the uniform (fast) components 
~La and Lib of the spins we arrive at the action of the 
nonlinear sigma model for the case K = (the detailed 
derivation of this model will appear elsewhere [12]). This 
action describes the low-energy, long-wavelength fluctu- 
ations of the staggered (slow) order parameter b (as 



usual we use /3 = 1/ksT and f = J dr J d 2 x) 
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The first line describes the low-energy, long-wavelength 
fluctuations of the two order parameters 11^,113, in- 
dependently. The second line contains their coupling, 
through 7, within the AF/Ising broken phase, where the 
term yia ■ is allowed (this model was derived for the 
(7r, 0) magnetic configuration, which breaks explicitly the 
Ising symmetry). For the Ising symmetric part of the 
phase diagram, such term would have been absent, but 
integration over fluctuations would give rise, instead, to 
a term like (n^ • n^) 2 , which is Ising invariant [l5| . The 
third line contains a dynamical term that describes the 
Bloch precession of the staggered moments in a given 
Neel state (say A) around the moments of the other (say 



B), and it resembles the coupling to a magnetic field. 
Finally, the fourth line contains interacting dynamical 
terms that also arise from the integration over La,b- 

All couplings in ([2]) are expressed in terms of the orig- 
inal parameters. The spin stiffness is p s = 2J2S 2 a 2 ~~ d , 
the spin- wave velocity equals Co = 2\f2Say/4J% — J 2 /h, 
we will also have a contribution from c\ = cq(J2/Ji), 
while the couplings between n^i and ns are given by 

7 = (l + 43*), = 4#7F> and ^ 

b = ^gr 1 (4j-i_ jz\ ■ The bare (unrenormalized) values of 
the parameters described above are, as expected, higher 
than the measured values. For example, for a = 5.695 A, 
Ji = 20 mcV,and J 2 = 40 meV, we find hc = 1.2 eV A, 
while the typical values are actually around 0.25 eV A. 
The theory described by Eq. ((2j) is highly interacting and 
a full renormalization procedure is required to reproduce 
the actual values of all c oup lings and the smallness of 
the Fe magnetic moment [l2j . For the purposes of com- 
parison with experiments in this letter, we shall use the 
already established values for some of these constants. 

To find the magnetic excitation spectrum we look at 
the poles of the staggered spectral function, A(k,u>). 
Within the Green's function formalism, these are ob- 
tained from the imaginary part of the retarded Green's 
function, G re t(k, w), for transverse staggered fluctuations 
,4(k,w) = -(l/7r)lina5_j.o^m[G ret (k,ia;„ -> ui + iS)], 
where u) n = 2Tm/ ft are the Matsubara frequencies. The 
first set of poles corresponds to 



huj(k) = Hcswyk 2 +jk x k 
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This is the usual acoustic 



where — . . 

c sw \ c o Ac \ 
spin- wave branch also found in linear spin wave theory 

[18j], which is gapless at the zone center, hu>(0) — 0, in 
agreement with Goldstone's theorem. Notice that the 
observable spin-wave velocity csw increases by a factor 
•y/2 + 7 as one moves from the k y — (or k x = 0) line 
to the k x — k y direction (along the spin stripes), show- 
ing that the spin- wave dispersion has a planar directional 
anisotropy, see Fig.[I] a consequence of the collinear char- 
acter of the antiferromagnetic order. The second set of 
poles corresponds to 



M7(k) = hcopJk 2 - jk x k y + 77, 



(4) 



where = [ \ + Here, instead, we find an opti- 

cal gap, Ml(0) = A op ^ 0, at the antiferromagnetic zone 
center (1, 0, 1), given by 



(5) 



which corresponds to long lived optical excitations. No- 
tice now that the optical spin-wave velocity along the 
k x = k y line decreases with respect to the k y = (or 
k x = 0) directions, see Fig. [TJ exactly the opposite from 
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Acoustic and Optical Branches (K-0) 
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AcusLic and Optical Gaps (K=0) 
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FIG. 1: (Color online): Acoustic (lower, black) and Optical 
(upper, red) spin-wave branches for the isotropic case, K = 0, 
as a function of |fc|, for two different directions, k y = (solid 
lines) and k x = k y (dashed lines). Inset: Acoustic and Optical 
gaps as a function of J1/2J2. Notice that while the acoustic 
gap is always zero, as imposed by Goldstone's theorem, the 
optical gap hardens as Ji increases, diverging for Ji = 2J2. 



the case for the acoustic branch. The optical gap hard- 
ens as J\ increases, see inset in Fig. [IJ and eventually 
diverges at J\ = 2J2, when the classical configuration 
changes and the action @ has to be modified. 

After including a single ion anisotropy (SIA) term 
at the Fe ions, K 7^ 0, we have extended the deriva- 
tion of the nonlinear si gma model accordingly (details 
will be given elsewhere 12]) and we have recalculated 
the magnetic excitation spectrum for K 7^ [12]. The 
rather lengthy expressions for the dispersions of the two 
branches for arbitrary K reduce, in the limit K <C J2, to 



fiw(k) = hcsw\/k 2 + jk x k y + A 2 SIA , 



hfl(k) = h Cop-Jit 2 - jk x k y + A 2 SIA + 77. 



(G) 
(7) 



As expected, the immediate effect of a SIA term is to 
produce a gap for all branches in the spectrum 
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which only exists for high spin systems, being absent 
when S = 1/2, like in the case of cuprates. 

We can now analyze recent inelastic neutron scattering 
results for spin excitations in SrFe2As2 (H)]. The spin- 
wave dispersion shown in Fig. [2j for scans along (H, 0, 1) 
at T = 160 K, can be well described by the lower branch 
of the spectrum, Huj(h), along the k x = k y line (we use 
a coordinate system rotated by 45° with respect to the 
original unit cell). This branch is gapped solely by the 
SIA and for T = 7 K the gap can be seen at 6.5 meV, 
see the inset a) of Fig. [21 However, as discussed in the 
introduction, the broadened integrated intensity of the 
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FIG. 2: (Color online): Spin- wave dispersion for Q = (H, 0, 1) 
and at T = 160 K. Experimental data from Ref. [131 ] . Inset 
a): integrated intensity at T = 7 K, as a function of energy, 
exhibiting a two-step profile consistent with the existence of 
two gapped branches in the magnetic spectrum, Eqs. (JSJ) and 
0. Inset b): Temperature dependence of the o ptic al gap, 
according to analysis of the data reported in Ref. [13l|. 
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Dirac Coties in Pnictidcs 
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FIG. 3: (Color online): a) (top left) Unfolded two-band Fermi 
surface; b) (top right) Folded four-band Fermi surface, result- 
ing from the staggered buckling of As ions; c) (bottom left) 
Dirac cones along the T — X direction, after folding of the 
Brillouin zone due to collinear AF order; d) (bottom right) 
Small hole pockets originating from the top of a Dirac cone. 



unpolarized neutrons, observed for either SrFe2As2 [13| 
and BaFe2As2 [14||, makes room for the existence of two 
or more gaps in the magnetic spectrum. In fact, the in- 
tegrated intensity shown in the inset a) of Fig. [5] shows a 
two-step profile, indicating the presence a second magnon 
gap at 10 meV. The temperature dependence of the op- 
tical gap, following a similar analysis of the data in |13| 
for T = 7, 80, 160 K, is shown in the inset b) of Fig. [1 

Next we analyze the quasiparticle excitation spectrum. 
The introduction of itinerant carriers into the problem 
presents us with an important question: if and how the 
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collinear AF order modify the quasiparticle spectrum. To 
investigate this problem we adopt a minimal two-band 
model for the quasiparticles, which can be described by 
the tight binding Hamiltonian 

K= E ttfKjid^ + h.c, (9) 

where (df^y creates an electron/hole at the site i, with 
orbital character a — d xz ,d yz and spin projection a. 
The values of the four hopping amplitudes, two for direct 
d xz ,yz-d xz , yz hopping and two for crossed d xz , yz -d yz _ xz , 
and d yz ^ xz — d XZt y Z , with different signs, and are the same 
as used in fill ] . Finally, the coupling between iron itiner- 
ant carriers to the iron local moments is done via flpj | 

H Hund = - J -f Y,Ka) f r(d?J ■ S it (10) 

i,a 

where f are the Pauli matrices and Jjj is the ferromag- 
netic Hund's exchange coupling constant. 

In order to calculate the quasiparticle spectrum we 
make use of a semiclassical approximation. We split the 
sum over lattice sites i in ()10|) into sums over sublat- 
tices ia and is, and we then replace Sa.b by their static 
equilibrium configuration within the colinear AF ordered 
state. The total quasiparticle Hamiltonian H qp + Huund 
is then quadratic in the (df CT ) T , (df a ) operators and can 
be easily diagonalized. The eigenvalues are plotted in 
Fig. [3] In the high temperature magnetically disor- 
dered phase the sublattice magnetization is zero and the 
dispersion is the one of the tight binding Hamiltonian 
H qp alone, where quasiparticles are characterized by a 
multi-band Fermi surface, with a rather large density 
of states, Figs. [3^) and b). Upon collinear AF order- 
ing, however, the doubling of the magnetic unit cell, and 
the consequent Brillouin zone folding, produces slightly 
anisotropic, hole-like Dirac cones at very specific loca- 
tions, (±0.257r, 0) and (±0.757r, 0), related to the collinear 
structure of the AF order, see Figs.[3t) and d), thus con- 
tributing a very small density of states. These two fea- 
tures are consistent with both recent ARPES results, see 
Fig. 1 of Ref. [19], and optical spectroscopy @, which 
reveal a strong reduction of the density of states upon 
AF ordering, for T < Tsdw (where T$dw is the spin- 
density- wave transition temperature). The existence of 
Dirac cones in iron pnictides had already been discussed 
in the literature [20( by using an itinerant-only model. 
Here, we have demonstrated that, also within a model 
where local moments interact with a bath of quasiparti- 
cles, such cones appear at the precise positions observed 
in experiments. Finally, let us remark that, from the 
point of view of the magnetic spectrum, the bath of car- 
riers simply provide a width (finite magnon lifetime) to 



the magnetic spectral lines, without otherwise modifying 
the pole structure obtained in this work. 

We have seen that the collinear AF ground state of 
the spin- Hamiltonian ([T]), containing two coupled, inter- 
penetrating Neel states, gives rise to a second, optical 
spin-wave branch in the magnetic excitation spectrum, 
that can be identified as an additional step in the in- 
tegrated neutron intensity in SeFe 2 As2 It remains 
to determine the specific polarization of each of the two 
spin-wave branches, through polarized neutron scatter- 
ing, and to look for the effects of such additional optical 
spin- wave branch in the thermodynamic properties of the 
SrFe2As2 system, such as the magnetic contribution to 
the specific heat [12j. We have also seen that, when cou- 
pled to a bath of itinerant quasiparticles, the collinear 
AF order produces hole-like, slightly anisotropic Dirac 
cones as recently observed by ARPES in Ref. [19j. 
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